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Abstract 

In the present paper, we have studied a discrete version of the WLC model, which incorporates 
the spatial constraints imposed by the magnetic tweezer, used in current micro-manipulation exper- 
iments. These obstruction effects are relevant for "short" molecules, involving about two thousand 

PQ 

base pairs or less. Two elements of the device have to be considered: first, the fixed plastic slab 

• 1— 1 \ 

on which is stuck one molecule end, second, a magnetic bead which is used to pull (or twist) the 
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attached molecule free end. We have developed quantitative arguments showing that the bead 
surface can be replaced by its tangent plane at the anchoring point, when it is close to the bead 
south pole relative to the pulling direction. We are, then, led to a confinement model involving two 



> 
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. repulsive plates: first, the fixed anchoring plate, second, a fluctuating plate, simulating the bead, 

O , 

( in thermal equilibrium with the attached molecule and the ambient fluid. The bead obstruction 



effect reduces to a slight upper shift of the elongation, about four times smaller and with the 



. same sign as the effect induced by the anchoring plate. This result, which may contradict naive 

expectations, has been qualitatively confirmed within the soluble "Gaussian" model for flexible 

> 

polymers. A study of the molecule elongation versus the contour length L exhibits a significant 
non-extensive behavior. Although the curve for "short" molecules is well fitted by a straight line, 
with its slope very close to the prediction of the standard WLC model, it does not pass through the 
origin, due the presence of an offset term independent of L. This leads to a 15% upward shift of the 
elongation for a 2 kbp molecule. Finally, the need for thorough analysis of the spatial constraints 
in super-coiled dsDNA elasticity measurements is illustrated by "hat" curves, giving the elongation 
versus the torque. 

PACS numbers: 87.15.By, 61.41+e 
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I. INTRODUCTION 



In the last decade, single particle biophysics has developed into a very active field of 



research 



5[. In particular, micro-manipulation experiments are now recognized 
as valuable tools to observe, in real time, a single double-strand DNA (dsDNA) molecule 
interacting with the proteins involved in the cell duplication processes. The basic principle 
is to look for sudden variations of the stretched dsDNA molecule elongation, which occur 
when the biochemical reaction is taking place. (For two recent reviews see the references 

In recent experiments [6|,|7|], there is a tendency to use relatively short segments with 2000 
base pairs (2 kbp), corresponding to L ~ 680 nm; this number is to be compared with the 
persistence length of the dsDNA molecule A ~ 50 nm. This implies that finite size effects 
may be of some importance, specially the spatial obstruction caused by magnetic tweezers. 

The simplest way to implement spatial constraints is to introduce in the dsDNA elas- 
tic energy density a one-monomer potential V(r(s)), where r(s) is the coordinate of the 
monomer, running along the chain of arc-length s. The "worm-like-chain" (WLC) model 
@, S 0] describes rather accurately dsDNA elongation experiments. In its usual formula- 
tion, the sole dynamical variable is the running tangent vector t(s) = 4-r(s) and in that 
case, one has to write r(s) = J S t(s')ds'. The potential energy to be added to E WLC takes 
then an ugly non local form J L ds V (J S t(s')ds'). This difficulty can be solved by formulating 
the model in such a way that t(s) and r(s) behave as independent dynamical variables. 

Numerous authors addressed this problem within the continuous version of the WLC 
model. The statistical properties of the molecular chain are obtained by solving a Quantum 
Mechanics problem, involving an imaginary time — i s. Using various arguments, they found 
that the Hamiltonian, allowing for spatial constraints, is obtained by adding two extra terms 
to the standard WLC Hamiltonian H WLC = — tj ~~ * ' f • (f ^ s the stretching force given 
in thermal units.) The first is the so-called "ballistic" term V ■ t and the second is the 



potential V(r) [lj, 
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24j . This model has been applied to 



various problems in semi-flexible polymer physics: the flow of semi-flexible polymers through 
cylindrical pores [3], the unbinding transition between semi- flexible polymers attracted 



by directe d p olymers [19] . the symmetric interface between two immiscible semi-flexible 
polymers 24( and probably others... 
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In reference jl3j , we have taken a different approach by remaining within a discrete version 
of the WLC model, which has to be introduced anyhow, if one wants to write down explicitly 
the functional integral giving the Boltzmann partition function. By using a simple trick, 
we were able to write the partition function as a multiple integral, where the coordinates 
r n and the tangent vectors t n of the discrete molecular chain are treated as independent 
integration variables. Using the transfer matrix formalism, it is then possible to write down 
a recurrence relation between adjacent intermediate partition functions Z n (r n ,t n ) relative 
to chains having a crystallographic length smaller than the actual one: 

Z n+ i(r n+ i,t n+ i) = exp(-6V(r n+ i)) J d 2 tt(t n )T WLC {t n+1 \ t n ) Z n (r n+1 - bt n+1 , t„) , (1) 

where TwLc{^n+i | t n ) is the transfer matrix relative to the unconstrained WLC model. All 
the explicit computations performed in reference [3] and in the present paper are based 
upon the above iterative construction, which has a suggestive interpretation in terms of a 
Markov random walk model in three dimensions. 

The confined dsDNA configurations studied in reference [3] were not fully realistic, since 
they do not account properly for the spatial obstructions occurring in magnetic tweezers. 
To get a feeling about the orders of magnitude involved in "short" molecules, say with 



L < 10 A, let us quote the results obtained in ref. [13j for the relative elongation upward 
shifts induced by the spatial obstruction of the anchoring plate. ( By anchoring plate, we 
mean the plastic slab upon which is stuck the initial molecular strand with the help of a 
"biological glue".) For a typical stretching force F ~ 0.3 pN, the upward shift is given by 
1.6 A/L, which amounts to 12 % for a 2 kbp molecule. In a magnetic tweezer, the free end 
of the dsDNA molecule is attached to a magnetic spherical bead, having a diameter of about 
one micron. In view of the above result, the bead obstruction effects are certainly worth 
investigating. (To get a very schematic view of the various constrained and unconstrained 
situations to be studied in the present work, see the upper graph of Fig [5] .) Such a study 
will be particularly relevant for dsDNA molecules with a few kbp, when they are stretched 
by forces within the range 0.1 to 0.5 pN. 

A theoretical analysis of the spatial obstructions in a magnetic tweezer is a difficult 
task, if one wants to treat it as a full three-dimension space problem, within the WLC 
model. The difficulties are about the same in the approaches based upon the solving of 
a Schrodinger-like equation or the transfer matrix iteration technique. The Monte-Carlo 
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method, which could perhaps be a viable approach, will not be considered in this paper 
because of lack of competence of the author. We are going to temper the above pessimistic 
views, by showing that under well defined conditions the nucleotides do not really feel the 
curvature of the bead. More precisely, the bead surface can be reasonably approximated 
by its tangent plane at the molecule free-end anchoring point, assumed to be the lowest 
point of the bead with respect to the force direction. It follows, then, that under realistic 
experimental conditions, the magnetic-tweezer obstruction can be simulated by two repulsive 
plates, normal to the stretching forces. There is, however, an important difference with the 
two-fixed-plate problem we have studied previously [13]. The initial molecular segment is 
still anchored to a fixed plate but the free end is now attached to a plate which is no longer 
fixed, but in thermal equilibrium with the dsDNA molecule and the ambient fluid. We are 
back to a one-dimension space iteration problem within transfer matrix technique. It is 
definitely more difficult than the two-fixed-plate problem but still manageable. 

Volume-exclusion effects in tethered-molecule experiments have been studied recently by 
Monte-Carlo techniques 25| . The authors have considered the situation where no stretching 
force are applied upon the bead. So, their significant work cannot be compared with the 
present one, since we are dealing with stretched molecules, having a relative elongation 
larger than 0.65. A finite-size effect analysis appears also in connection with the entropic 



261 ] . In particular, the authors deal 



elasticity of DNA molecule, having a permanent kink 
with the boundary conditions to be satisfied by the tangent vectors at the two ends of the 
molecular chain. In the present work, we have concentrated on the spatial confinement, 
imposed upon the internal monomers, by the repulsive surfaces holding the two free ends. 
It is easily seen that these constraints insure that the initial and terminal tangents vectors 



do satisfy automatically the "half-constrained" boundary conditions of ref. 



26]. 



II. A SIMPLIFIED MODEL TO DESCRIBE THE MAGNETIC BEAD SPATIAL 
OBSTRUCTION. 

In this section, we would like to develop arguments to justify the replacement of the 
bead surface by its tangent plane at the anchoring point, assumed to be located at the bead 
south pole with respect to the force direction. The discussion will be performed within a 
discrete version of the WLC model. The molecular chain is represented by N elementary 
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links, involving point-like "effective" monomers, separated by a length b, much smaller than 
the persistence length A. The effective monomer number N is related to the contour length 
L by the relation N = L/b. A microscopic state of our model is then defined by the set of 
2N vector variables: {(r^tx) ....(r n , t„), ...(rjy, tjv)} with 1 < n < N, where r n and t n , are 
respectively the monomer coordinate and the unit tangent vector such that bt n — r n — r n _i. 
In this paper the z axis is parallel to the direction of the stretching force F and has its 
origin at the fixed-plate anchoring point. One must stress that the variables and r n with 
1 < n < N — 1 have to be treated on different footing: is the coordinate of the terminal 
monomer but it gives also the position of the bead. 

We shall, first, discuss the transverse fluctuations of the terminal monomer, (x 2 N ); they 
are given in the thermodynamic limit L/i > I by a well known formula 5] : 

(x > N ) = (z N )(k B T)/F = ^A, (2) 

a 

where a is the dimensionless force parameter F A/iksT). This formula is used to calibrate 
the force in magnetic tweezer experiments and was first obtained by a simple thermodynamic 



argument Q]. In an unpublished note- 



Ill , we have given a statistical mechanics evaluation of 
(x 2 N ) within the WLC Model. We have recovered the formula in the limit L/A ^> 1 with 
corrections of the order of A/L, which cannot be totally ignored for the "short" molecules 
considered in the present paper. For this reason, we have given an updated version of this 
note in the appendix. The above value of (x 2 N ) will be used as a benchmark in our estimate of 
the transverse fluctuations of the internal monomers. To proceed it is convenient to introduce 
cylindrical coordinates for the effective monomer positions r n = (z n , rj_ n ) and the associated 
tangent vector t n = (cos # n , tj_ n ). Our purpose is now to estimate the thermal average 
of the square of the transverse distance between the n-monomer and terminal monomer: 
( ( r ±n ~~ t ±n) 2 )- More precisely we are going to establish the inequality: 

( (rj. B - r ±N ) 2 ) < ( v ± 2 N ) = 2( x 2 N ) = 2^ A. (3) 

a 

We have found convenient here to work within the WLC model in its simplest form, where 
only the tangent vectors t n = (cos 9 n , cos# n , sin# n cos <p n , sin 6 n sin0 n ,) appear explicitly. 
The n-monomer coordinate is then given by the sum: r n = b £)i=i tj. The discrete WLC 
model is best formulated in terms of the transfer matrix connecting two adjacent links: 

/ a 

Twlc(K+i, tn) oc exp (-— (t n+ i - t n ) 2 



oc exp fy(cos0 n+ i cos# n + sm9 n+ i sin9 n cos(0 n+ i - (/> n ))) , (4) 

where for simplicity we have omitted the stretching energy term. Our aim is to prove that 
the following difference is positive: 

i=N j=n 

A ± = (r 1 _ 2 N )-((r ±n -r ±N y) = (r ± l) + 2 £ £<t ±l -t ±J >. (5) 

i=n+l j=l 

Clearly, we have to show that the thermal average: (t_y • t±j) = (sin#j sin#j cos(0j — 4>j)) 
is positive. Within the usual definition of the cylindrical coordinates, sin 6i sin 6j takes only 
positive values while cos(0j — 4>j) can be negative as well as positive. To prove the positivity 
°f (t_u ' t_L,-}, it is then enough to perform the thermal average over the azimuthal angle 
fa, the longitudinal components cos#j being frozen. Assuming that fa is initially uniformly 
distributed, the overall system is invariant upon any global rotation around the z-axis, so 
that we can replace the set of the N variables fa by the N — 1 variables ipi — fa — fa-i with 
1 < i < N. Using the formula (jlj) one gets immediately the probability distribution of fa : 

Pi{i>i) = exp (—Xicosipi^ I (2-7rio(Ai)) , A» = sinftj sin^j > 0, (6) 

where, for convenience, we shall take — ir < ipi < tt. We have, now, everything we need to 
compute (cos(0j — fa))- We first note that fa — <f)j can be easily written as a sum of ipi 



angles: 



l=i 

5>i-i + fa-i - fa-2--fa+i ~fa= ^i- 



1=3+1 

The average (cos(0j — <f)j)) is then easily performed: 

l=i l=i 

(cos(fa - fa)) = Re( Yl (exp(#;)) = II ( cos ^z)> ( 7 ) 

l=j+l l=j+l 

/ — 7T 
dipiPi(ipi) cosipi. (8) 

By looking at the functional form of Pi{ipi) in eq. it is easily seen that Xi > implies 
(cos tpi) > and, as a consequence, the positivity of (cos(0j — 4>j))- Since the final averaging 
will preserve this positivity, we can conclude from eq. ([!]) that > 0, or in a more concrete 
way: 

<(r^-r^) 2 ><(r^) = ^A (9) 

a 

Let us consider, now, an internal effective monomer going upward and crossing the tangent 
plane at a transverse distance |(rj_ n — r± N )\ from the south pole. The maximum vertical 



distance 5 curv which it can travel before hitting the bead surface is given by a simple geo- 
metrical argument: 5 curv = (rj_ n — r± N ) 2 /(2 R) where R is the bead radius. Performing the 
thermal averaging we get the final inequality: 

(^™)<^ = ^A (10) 

The internal monomers are subjected, in the vicinity of the bead, to the chain tension 
force which tends to pull them towards the anchoring point. As a consequence, the above 
assumption of a vertical path leads to an overestimate of S curv , so relaxing this constraint 
can only but strengthen the inequality of eq. fllOl) . 

In order to quantify the internal monomer ability to detect the curvature of the bead, 
we are going to compare 5 curv to other experimental lengths. Let us begin by the average 
elongation: (z(N)). We readily obtain the inequality: 

(Scurv) /(Z(N)) <4~ = L25 X 10 ~ 2 AR(/4 (11) 
HO. 

where we have taken A = 50 nm and a = 4, which corresponds to a stretching force 
F = 0.31 pN. R(fi) stands for the bead radius given in micron. A more significant com- 



parison involves the mean-square longitudinal fluctuation Az^ — J (z 2 N — (zn) 2 )- We have 
computed the ratio Az^/ (zjf) using a version of the WLC model incorporating the confining 



effect of the fixed plate holding the initial strand [13j. Taking a = 4 and L = 12 A, we have 
obtained: A%/ (zn) = 0.083. Keeping the same value of L, we have derived upper bounds 
of the ratio (5 curv ) / A zn, for a = 2 ,4, 5 respectively: 

(Scurv) / A z N < V^-TT- = (°' 17 > °' 15 ' 0.13)/i?(/x). (12) 
A^tv Rot 

We conclude that the internal monomer mean free path above the tangent plane (S curv ) is less 
than one sixth of the mean-square longitudinal fluctuation of the bead A z^ when L = 12 A 
and a > 2. This suggests that, under such conditions, the effective internal monomers are 
not really able to detect the bead curvature and gives a rather strong justification for the 
replacement of the bead surface by its tangent plane at the anchoring point. 

However, one must keep in mind that the above conclusion hinges upon the simplifying 
assumption that the end of the DNA molecule is stuck at the south pole of the bead. This 
condition will not be satisfied, unless some selection is performed among the beads. This can 
be done in practice by slowly rotating the magnetic tweezer around its axis. The magnetic 




FIG. 1: A schematic picture showing the basic physical parameters governing the motion of an 
internal monomer (small black dot) in the vicinity of the bead surface tangent plane at the terminal 
monomer ( big black dot) anchoring point at the south pole. The figure represents a section of the 
bead by a plane containing its center and parallel to the z-axis. The magnitude of the fluctuations 
are shown in the case of a force F = 0.31 pN and a contour length L = 12 A. As indicated in 
the text, the yellow vertical band gives an upper limit to the allowed transverse motion of an 
internal monomer with respect to the terminal one. The red region represents the vertical section 
of the excluded volume coming from the replacement of the bead surface (the blue circle arc) by its 
tangent plane. The emerald horizontal band gives the amplitude of the longitudinal fluctuations 
of the bead. A comparison with the vertical width of the red region suggests that the internal 
monomers will not be able to "feel" the bead curvature. The big blue dot represents the case 
of an anchoring point lying away from the south pole. The plane simulating the bead is now 
the horizontal plane z = (its section appears as a blue dashed line). The excluded volume is 
subjected to a positive variation on the right-hand side and a negative one on the other side; so 
there is clearly a compensation effect. Furthermore, the selection procedure described in the text is 
affected by the transverse fluctuations of the bead. On this figure, we have assumed for simplicity 
that the selection has been performed with F = 0.31 pN. Taking instead F = 2.7pN would have 
reduced the selection angle 9u m by a factor 3. 



bead behaves as a compass and follows the rotating magnetic field. If the experiment is 
performed at a high enough force F se i ec , the bead rotates around the vertical axis passing 
through the anchoring point lying away from the south pole ( the big blue dot appearing 
on Fig{T]). As a consequence, the center of the bead describes a circle of radius RsmO an , 
where 6 an is the "latitude" of the anchoring point with respect to the bead south pole (i.e 
the angular distance between the black point and the blue point on Fig{T]). There is clearly 
a limitation in this bead selection procedure, coming from the transverse fluctuations of the 
bead. Only the beads satisfying the inequality sin9 an > sin6^ m = J (r±%) (F/ F se i ec ) j R 
can be eliminated in practice. (We have made explicit the fact that the actual elongation 
experiment is done at a force F different from that used in the selection procedure). Let us 
consider now a bead belonging to the selection. Replacing the tangent plane by the horizontal 
plane z = induces a variation (positive or negative ) of the vertical free path of the internal 
monomers of amplitude smaller than sin6^ m J (r_|_^) = J F j F sdec (r±j^)/R. Taking into 
account the compensation apparent on FigJT] and the eventual reduction factor \j Fj F se \ ec 
the average height of the excluded volume can easily be made one order of magnitude smaller 
than the longitudinal fluctuation of the bead Azn- Our model is then expected to be valid 
for the beads selected according to the above criterion. 



III. A SOLUTION OF THE MODEL AND ITS PHYSICAL INTERPRETATION. 



Computation Procedure. 



13fj to the 



In this section, we are going to extend the matrix transfer method of reference 
study of the statistical properties of a dsDNA molecule confined within a magnetic tweezer. 
The two ends of the molecule are respectively attached to a fixed plate and to a magnetic 
bead immersed in a liquid which simulates the cellular medium. The bead is subjected to an 
external force normal to the anchoring plate. This pulling force is balanced by the tension of 
the stretched molecule. As shown in the above section, the identification of the bead surface 
to its tangent plane at the terminal monomer position is a fairly good approximation if the 
anchoring point is close to the south pole. 

Our confinement configuration is then defined by two parallel repulsive plates, the first - 
the anchoring plate - is fixed and the second - simulating the bead - is in thermal equilibrium 
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with the attached molecule and the ambient fluid. With proper initial conditions, the 
partition function of the system molecule-plus-bead is invariant, first, under rotations around 
the stretching force direction and, second, under translations parallel to the plates. We are 
going to use the same discrete version of the WLC model as in the previous section. A 
microscopic state of our model is then defined by the set of 2N longitudinal variables: 
{{z±, 9x) ....(z„, 9 n ) ...{z N) 9 N )} with 1 < n < N. 

We proceed in two steps. We, first, assume that the terminal monomer coordinate z^ 
has a fixed value, taken among a finite set chosen to be representative of the actual physical 
spectrum. The internal monomer coordinates z n , with 1 < n < N — 1, are truly stochastic 
variables, associated with the partition function Z n (z n , 9 n \z]y). For the moment, z^ is treated 
as an external physical parameter. When 1 < n < N — 1, the internal partition functions 
Z n {z n , 9 u \zn) obey a recurrence relation which is just a rewriting of a formula given in 
Section 2.1 of reference [3]. The rule is very simple: take the iteration equation for the 
fixed-plate confinement problem and identify zn with the distance between the two repulsive 
plates L : 

Z n (z n ,9 n \z N ) = exp(-bV(z n ,z N )) I ' d(cos9 n ^ 1 )T WLC (9 n ,9 n ^ 1 , f) x 

Jo 

Zn-\{z n ~ b cos 9 n ,9 n - 1 \z N ). (13) 

The potential V(z n , z^), is written as the sum of two terms: 

V(z N , z n ) = V p i(z n ) + V p i(z N - z n ), (14) 

where the first one is associated with the fixed anchoring plate and the second with the 
fluctuating plate simulating the magnetic bead surface. The repulsive plate potential V p i{z) 
is given in terms of the rounded-off step function 131 ]: 

&(z,Az) = ± + ±erf(z/Az) } (15) 

where erf (x) is the "error" function : -^=/o ex P( — t 2 ) dt and Az a smoothing length assumed 
to be ~ A. For the sake of simplicity, it is convenient to give directly the associated 
Boltzmann factor: 

ex V (-bV pl (z)) = e(z,Az). (16) 

The conditional probability distribution T WLC {9 ni 6* n _i, /) is obtained by performing an az- 
imuthal average of the transfer matrix given by equation (jlj): 

A b f 

TwLc(9i,9 2 ,f) = exp-{— (1 - cos#i cos# 2 ) + — (cos#i + cos0 2 )} 
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A sin 6>i sin 9 2 

x^o( 1 ), (17) 

where / is related to the stretching force by F = f k^T. We compute the partition function 
relative to the last internal monomer, for a fixed terminal monomer, Z^^x{z^^x,6 n \zjs[), by 
running, up to n = N — 1, the iteration process defined by equation ( |T3l) . To perform the 
relevant recurrence process we have used Mathematica codes where analytical and numerical 
computations are intertwined. For more details about our procedure, see the SectionV of 
reference 29| . 



The second step ( a short one !) is to compute the partition function Z N (z^, 9^) relative 
to the terminal monomer. It is easily obtained from the formula: 

Z N (z N ,9 N ) = exp(-bV p i(z N )) f d(cosB N ^ 1 )T WLC (6 N ,6 N -i, f) x 

JO 

Z N -i{z N -b cos9 N ,9 N „ 1 \z N ). (18) 

Finally, the probability distribution relative to the terminal monomer longitudinal coordinate 
- it is also the longitudinal distance of the bead south pole from the fixed plate - is given as 
follows: 



1 f 1 

Pn(zn) = jj= J^d(cos6 N ) Z N (z N ,9 N ), 

j\f = f L d(z N ) ( l d{cos9 N )Z N {z N ,9 N ). (19) 
Jo Jo 



B. Result and Physical Interpretation. 

The probability distribution Pn(zn) is displayed upon the left graph of Figj2]as a contin- 
uous blue curve. It appears together with green and red curves, corresponding, respectively, 
to the two situations: absence of constraints and presence of a fixed anchoring plate. Two 
features are conspicuous: 

• The red and blue curves are very close, while the green and red ones are much further 
apart, showing that the obstruction effect of the bead is much smaller than that coming 
from the fixed anchoring plate. 

• Perhaps even more surprising, the bead obstruction effect reduces to a slight push 
upwards given to the terminal monomer, while one might have expected naively the 
reverse. 
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In the two-fixed-plate configuration studied in our previous work [13|], the confinement 
effects upon Pn(zn) were very spectacular ( see Figure 4. of this reference). It must be 
stressed that the physics involved was very different from that prevailing in the present 
paper. Indeed, the terminal monomer, supposed to be attached to a nano- magnet, was free 
to move between two diamagnetic fixed plates. When the elongation in absence of plates 
was larger than the two-plate distance, the probability distribution Pn{zn) was flattened 
against the upper plate by the stretching force. No such phenomenon appears here since the 
terminal monomer is stuck upon the plate simulating the magnetic bead. Only the internal 
monomers feel the repulsion of the magnetic bead and it is through their intermediary that 
the terminal monomer is affected by the spatial obstruction of the bead. 

A plausible mechanism goes as follows: the internal monomers are pulled upwards by the 
string tension force. When they collide with the bead, they are expected to give a small 
upward push to the bead surface. In other words, they are exerting an upward pressure on 
the bead, which is transmitted to the terminal monomer, leading to a small increase of the 
stretching force. 

The positivity of the variation of the elongation S(zn) can also be obtained by a simple 
thermodynamic argument, valid in the limit S(z^) / (zn) "C 1. The internal energy of the 
bead-plus-DNA system receives a positive contribution 5U coming from the repulsive inter- 
actions: Y.n=i Vpi( z N — z n) > 0. Writing that the variation of the free energy T = U—(z N ) F 
vanishes at equilibrium, one gets: 6(zn) — SU/F > 0. 

In order to get a confirmation of the above picture, we have repeated our computations 
within the soluble "Gaussian Model" , often used to describe "flexible" polymers. It involves 
a chain of point-like monomers connected by harmonic springs. The continuous limit of the 
chain is a string described by the elastic energy linear density: S 9aus (s) = ± (r(s)) 2 + V(r), 
where r(s) is the derivative of the monomer coordinate r(s) with respect to the string arc- 
length s and a -1 is proportional to the rigidity of the spring connecting nearest-neighbor 
monomers. The value chosen for the parameter a guarantees that the relative elongation af, 
in absence of spatial constraints, is the same as in the WLC model. The monomer number n 
is related to the arc-length s by s = nh with b/a ^ 1. If we identify V(r) with the potential 
V(z n ,z^) introduced above, the internal monomer probability distribution, P^ n *(x„, y n , z n ), 
factorizes into three independent distributions relative to each component. In the continuous 
limit, the statistical properties of the internal monomers are easily obtained by exploiting 
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FIG. 2: Probability distributions of the terminal coordinate zn, relative to different spatial ob- 
struction conditions, as predicted by two models. The continuous curves appearing on the left-hand 
lower graph have been obtained within the WLC model, while the dashed ones on the right-hand 
graph are relative to the Gaussian model. As indicated in the picture given in the upper graph, 
the green color corresponds to a molecular chain free of any external spatial constraints. The red 
and blue colors are associated respectively with two kinds of spatial obstruction: restriction to 
the upper half space by an anchoring fixed plate and confinement between the fixed plate and the 
fluctuating bead, holding the terminal monomer. 

the analogy with a QM problem involving the following simple Hamiltonian: 



The Hamiltonian relative to the terminal monomer H term is obtained by performing the 




a d 2 



(20) 



replacement : z, 



n 




The probability distribution for the ter- 
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minal monomer involving a fixed anchoring plate together with a pulling one, in thermal 
equilibrium with the "Gaussian" polymer chain, reads then as follows: 

Pn(zn) = J dz N -! exp(fz N )(z N \exp(—b # term ) \z N _ 1 )(z N _ 1 \ exp(— 6 (N — 1) H int )\z ). 

(21) 

The results obtained within the "Gaussian Model" are displayed on the right-hand graph of 
Figj2j The obstruction effects are qualitatively similar to those obtained within the WLC 
model but the chain-elongation increase induced by "the bead" 5 (z^) is much larger. This 
amplification can be understood by noting that the "Gaussian flexible" monomers, being 
allowed to move much more freely, have a larger collision rate with "the bead" . 

We have studied the magnetic tweezer spatial constraint effects on the dsDNA elongation 
within the stretching force range defined by 2 < a < 5, in the case of a contour length 
L = 12 A. The elongation curves are displayed in Fig [3] for three configurations : no 
space constraints ( green line), an anchoring plate barrier ( red line) and an anchoring plate 
barrier together with the magnetic bead ( blue line ). As discussed above, spatial obstruction 
effects from the bead, in thermal equilibrium, lead to a small increase - about few %- of the 
elongation. It is to be compared with the four times larger effect induced by the anchoring 
fixed plate, where the repulsive character of the fixed barrier is clearly exhibited. 

IV. THE NON-EXTENSIVE BEHAVIOR OF THE DNA ELONGATION IN 
PRESENCE OF SPATIAL CONSTRAINTS. 

In this section we shall study, for fixed forces, the variation of the elongation (z(N)) 
with respect to the contour length L. This analysis has been performed previously in ref. 
{"Lsl ]. In the case of fixed plates, it was found that (z(N)) is no longer an extensive variable 
with respect to L. The effects were very spectacular in the two plates configuration when 
the distance L between the two plates is smaller than the elongation in absence of spatial 
constraints. Conversely, when Lq > L, the molecule feels only the anchoring plate and the 
extensive behavior of (z(N)) is perturbed in a very simple way. For stretching forces such 
that a > 2 the derivative of elongation d/dL (z(N)) does not vary with L and stays very 
close to the constant (z(N)/L) WLC predicted, for a given force, by the WLC model in 
absence of spatial constraints. The only modification is the apparition in the elongation 
of an offset term independent of L. In more precise words, the elongation (z(N)) can be 

14 



< z >/L l= 12 A 




FIG. 3: This figure displays three curves, giving the "short" dsDNA relative elongation versus the 
stretching force, corresponding to different spatial constraints. Going from top to bottom, the blue 
curve is the prediction of the WLC model, incorporating both the bead and the anchoring plate 
spatial obstruction. The red curve accounts for the sole effect of the anchoring plate. The green 
curve, lying significantly below, has been computed in absence of any spatial constraints. These 
curves emphasize the dominant role played by the anchoring plate. 

written as follows when L > 2 A: 

(z(N)) = L (z(N)/L) WLC (1 + e) + AA f s (a), (22) 

where |e| < 1CT 2 and the dimensionless offset A ^ s (a) is a slowly decreasing function of a, 
which takes values of the order of unity when 2 < a < 5. 

In the present paper, we have performed - within our model - the same analysis in presence 
of the magnetic bead in thermal equilibrium with the dsDNA. We have also found a linear 
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FIG. 4: dsDNA elongation versus contour length, within the interval 2 A < L < 12 A. The 
blue points have been obtained from a version of the WLC model, implementing the magnetic 
tweezer spatial obstructions. They are well fitted by a straight line (the blue line), having a slope 
very close to the prediction of the unconstrained WLC model, given by the green line. It is to 
be noted that the extrapolated fitted line does not pass through the origin, as it should if the 
corrected elongation were still an extensive quantity. The constant offset term, responsible for this 
non-extensive behavior, gives rise to a 15% upward shift of the elongation when L = 12 A. 

variation of the elongation with respect to L, similar to that given by equation (I22P ; the 
only difference is an increase of the offset function A Q f s (a) by a few tens of percents. As 
an illustration, we have plotted on Fig. H] the result of a linear fit (the blue line) involving 
the elongations (z(N)) relative to a given set of contour lengths. They were computed with 
the method described in section II for a = 4, when L = N b takes 10 equally spaced values 
within the interval 2 A < L < 12 A. The slope coming from the fit coincides, to better than 
1%, with the prediction of the unconstrained WLC model, which appears as a green line 
on Fig. HI Similar features hold true for the linear elongation fits performed on the results 
obtained with seven equally spaced values of a within the interval 2 < a < 5. Although we 
have verified the validity of the linear fit for values of L > 12A and a > 2 , it is certainly not 
valid for a C 2. We have, indeed, found in ref. [13| (section 2.2) that the elongation upward 
shift contains an extra term cx a/L when < a < 1, with a coefficient which exhibits a very 
steep decrease with a. 
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FIG. 5: Numerical values of the offset term A f s (a) = (1 — L /dL) (z(L))/A versus seven rep- 
resentative values of a. A second order polynomial involving a -1 provides a rather accurate 
interpolation within the interval 2 < a < 6. 

We have plotted on Figj5] the values of the offset function A f s (a) coming from the 
corresponding elongation linear fits. For the sake of convenience, we have performed a third 
order polynomial fit, using a -1 as variable: 



No particular physical significance is to be attributed to the functional form chosen for the fit 
other than the fact it gives a rather simple and accurate interpolation between the calculated 
values. The formulae (1221) and ( BSD , together with the accurate values for (z(N)/ L) WLC 
given in reference [h]], lead to a one % evaluation of the dsDNA elongation, corrected for 
the spatial constraints induced both by the anchoring plate and the magnetic bead. A word 
of caution: formula (|23|) is not to be trusted if it is used outside the range 2 < a < 5 and 
for L < 2 A; in particular, it does not give the correct limit for a ^> 1 since one expects, on 
physical ground, that lim^oo A f s (a) = 0. 



A o/a (a) = 1.6942 + 0.8052 a" 1 + 0.131213 a 



-2 



(23) 
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V. THE EMPIRICAL DETERMINATION OF THE PERSISTENCE LENGTH 
FOR SPATIALLY CONFINED "SHORT" DSDNA MOLECULES. 



One may consider that the modification of the dsDNA elongation induced by the magnetic 
tweezer obstruction effects is, after all, not so dramatic since for the "typical" case a = 
4 and A/L = 1/12, it amounts to an increase of 15%. However, if the elongation together 
with transverse fluctuation measurements are used to get a determination of the persistence 
length, the situation becomes much more serious. We would like to show that the raw data 
of (z(L)), without a subtraction of the offset correction A f s (a)A/ L, could lead to totally 
unreliable values of A. The basic principle behind this determination of A is to combine two 
formulae giving the stretching force F in terms of {z(L)) and (x 2 (L)): the first one is the 
Strick et al. |5] formula, which is discussed in details in the appendix; the second one is the 
relation between the reduced force parameter a and the relative elongation u = (z(L))/L. 
In the case of the unconstrained WLC model, it takes a scale invariant form: a = T{u), 
where J-{u) is a numerical function, which interpolates accurate numerical results [l2]. By 
a simple algebraic manipulation, one arrives at the relevant formula, valid only within the 
unconstrained WLC model: 

It cannot be applied bluntly to the case of "short" DNA molecules studied in the present 
paper. Let us call u raw = (z(L)) raw /L the relative elongation involving the non-subtracted 
elongation (z(L)) raw , given in our model by equation (1221) and A raw the corresponding per- 
sistence length, obtained by plugging u raw in equation (I2~4j) . Ignoring the effect of spatial 
constraints upon the transverse fluctuations (x 2 (L)), one arrives immediately to the ratio: 

■rt-raw/rt-— sp, \ * 

Performing the numerical evaluation for the "typical" case a = 4 and A/L = 1/12, one gets 
the rather spectacular result: A raw j A = 2.9. Such a large number cannot be be explained 
by the finite size correction to the Strick et al. formula, estimated in the appendix. It is 
coming from the fact that J-{u) exhibits a very sharp increase when u > 0.75, due to the 



presence of a pole singularity at u — 1. (See the foot note [12J.) If one subtracts from u raw 
the offset correction A Q f s (a)A/ L, the ratio gets back to a value very close to 1. 

If the non-extensive behavior of the elongation has the simple linear behavior described in 
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the previous section, the required subtraction may be achieved empirically. The procedure 
involves elongation measurements upon three (or more !) "short" molecules with different 
contour lengths, say Li = 5 A , L 2 = 10 A and L 3 = 15 A, subjected to the same stretching 
force. (This latter requirement may be difficult to satisfy with precision.) A linear fit to 
the data will lead to an empirical determination of A A f s (a), together with a verification 
of the near equality between the fitted slope and the prediction of the unconstrained WLC 
model. If the result of the latter test is positive, then one can proceed to the required 
subtraction from the elongation and plug the result in equation (124)) . One should get in this 
way a value of A close to that obtained for "long" molecules, say with L/A > 100. Finally, 
the comparison of the empirical value of A f s (a) with the prediction given in the previous 
section will provide a further significant test of our model for the spatial obstruction effects 
in a magnetic tweezer. 



VI. SUMMARY AND PERSPECTIVES 



In the present paper, we have proposed a theoretical model for the spatial-constraint 
corrections to the dsDNA elasticity measured with a magnetic tweezer. An evaluation of 
the obstruction effects of the fixed anchoring plate had been given previously 13j and here we 
have concentrated upon the magnetic bead which is attached to the molecule free end. The 
main and somewhat unexpected result of this work is that the magnetic bead obstruction 
effects give rise to a slight upward shift of the elongation, about four times smaller than the 
anchoring plate effect. 



A. Synopsis of the paper. 

• In section II, we have developed theoretical arguments to justify the replacement of 
the bead surface by its tangent plane at the anchoring point, assumed to be located 
close to the bead south pole, defined by the force direction. As a first step, we prove, 
within a discrete version of the WLC model, that the mean square transverse distance 
between an arbitrary internal "effective" monomer and the anchoring point is smaller 
than the transverse fluctuations (^ir) of the terminal monomer. This latter quantity is 
given by the Strick et al. formula |5j in terms of the force and the molecular elongation. 
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Then we proceed, by a simple geometrical argument, to the derivation of a lower bound 
for the internal monomer mean free path 5 curv above the tangent plane. For "short" 
molecules ( L ^ 10 A), stretched by a force F > 0.3 pN, 5 curv is about six times smaller 
than the root mean square longitudinal fluctuations of the terminal monomer, if the 
bead radius is larger than one micron. This result suggests strongly that the internal 
monomers do not really "feel" the bead curvature. 

In section III, we have given a transfer matrix solution of our confinement model, 
involving a discrete chain of iV effective monomers with the two extremities attached 
to a pair of parallel repulsive plates. The initial monomer is anchored upon a fixed 
plate, while the terminal one is stuck to a fluctuating plate, in thermal equilibrium 
with the attached chain and the ambient fluid. In the first step, the terminal monomer 
longitudinal coordinate has a fixed value zn, taken among a representative set. The 
partition function of the chain of N—l fluctuating internal monomers are then obtained 
by solving, for each fixed value of z^, a two-fixed-plate confining model, using the 
method of ref. [3]. In the second step, we obtain the partition function of the terminal 
monomer by applying to the last internal monomer the transfer matrix relative to the 
spatially unconstrained WLC model. The probability distribution of the terminal 
monomer distribution Pn(Zn) exhibits two remarkable features, when it is compared 
to the same distribution in absence of magnetic bead obstruction. First, the two 
curves, appearing in Fig. 1 as blue and red continuous lines, are very close. Second, 
the effect of the bead reduces to a slight upward push given to the terminal monomer. 
One may have expected naively exactly the reverse: a significant downward push 



13j], only 



from the "bead". However, in contrast with the two-fixed-plate problem 
the internal monomers feel directly the repulsion from the "bead" , while the terminal 
monomer stuck on the bead don't. So, one has to invoke an indirect effect involving 
the internal monomers. A plausible candidate is the upward pressure they exert on 
the bead, which is transmitted to the terminal monomer. As long as we are dealing 
with a small effect, a simple thermodynamic argument leads also to an upward shift 
of the elongation. We have obtained a qualitative confirmation of the whole picture 
within the soluble "Gaussian" model, used for flexible polymers. We have found similar 
effects, but about seven times larger in typical cases. This amplification reflects the 
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fact that the internal monomers of a flexible polymer, being allowed to move more 
freely, have larger colliding rates with the "bead" . 

• The section IV is devoted to an analysis of the non-extensive behavior of the DNA 
elongation induced by the magnetic tweezer confinement effects. We have, indeed, 
found that in DNA molecules, having a contour length L > 2 A, the elongation (z(N)) 
is no longer an extensive quantity, within the force range 2 < a < 5. The derivative of 
the elongation with respect to L, d/d L(z(N)) , stays very close to the constant value 
predicted by the unconstrained WLC model. The sole non-extensive effect is the 
apparition in (z(N)) of an offset term A Q f s (a) A independent of L. In other words, 
the elongation-versus-L curve is still a straight line but it does not go through the 
origin. The dimensionless offset term A / S (a) decreases slowly from 1.6 to 1.3 within 
the interval 2 < a < 5 and it is well represented by a second order polynomial in 
a . For the "typical" case a = 4 — > F = 0.3 pN and A/L = 1/12, the non-extensive 
correction to (z(N))/N amounts to 15%. 

• In the final section, we investigate the possible influence of the magnetic tweezer con- 
finement effects upon the determination of the persistence length A, from elongation 
and transverse fluctuations measurements. We consider the case of "short" molecules 
of about 2 kbp, stretched by a force F ~ 0.3 pN. Plugging in the basic formula a "raw" 
relative elongation, uncorrected for non-extensive effects, leads to an overestimate of 
A by a factor 3 with respect to "long" molecule values. We suggest an empirical way 
to perform the required subtraction from the measured elongation. 



B. Possible extension of the present work to super-coiled dsDNA molecules. 

In references 
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29j the WLC model has been generalized to a Rod Like Chain 
(RLC) Model, involving both bending and twisting rigidities. This makes possible the study 
of super-coiled dsDNA entropic elasticity below the denaturation threshold. 

One can readily modify the RLC model in order to incorporate spatial constraints. The 
recurrence relation for the partition function Z n (z n ,8 n , k) relative to a super-coiled DNA 
molecule, with a given torque V = ksTn acting upon its free end, is obtained by per- 
forming in the r.h.s. of the recurrence relation (fl3l) the replacement: 7wLc(8 n+ i,9 n , f) — > 
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FIG. 6: Two "hat" curves, giving the elongation of a super-coiled molecule versus the torque. The 
upper curve incorporates the obstruction effects of the anchoring plate, which are ignored in the 
lower one. 



1~HLc(@n+i,@n, f, — ^ 2 ), where T RLC is given explicitly in ref. 29]. The anchoring-plate bar- 
rier is expected to have significant effects upon the so-called "hat curves" , giving, for a fixed 
force, the relative elongation versus the super-coiling reduced parameter a. Let us take the 



"low" force case studied in ref. [28| , F ~ 0.1 pN, where the RLC model "hat curve" dips 
steeply into the negative z region when \a\ > 0.03. This effect is attributed to the cre- 
ation of plectonem structures which are allowed to wander in the z < half space, because 
of the vanishing of their stretching-potential energy. Therefore, we can expect important 
modifications of the torsion elasticity if the spatial constraints are incorporated in the RLC 
model. 

This effect is illustrated in Fig. 5, which displays the results of a preliminary computation 
giving the relative elongation (Z(L)) versus the torque T = nksT for a "short" DNA 
molecule ( L/A=12) and F = 0.33 pN. The solid line is a spline fit connecting the points 
obtained with the RLC model, taking into account the anchoring plate obstruction effects. 
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The dashed line is relative to the RLC model free of any spatial constraints. The dip of the 
latter curve into the z < half-space reflects the formation of plectonem-like configurations. 
For < k < 1, the plate barrier effect reduces to an upward shift of the elongation, similar 
to that found in absence of super-coiling. For < k < 1 the anchoring plate confines the 
molecular chain in the z > half space. 



APPENDIX: CALIBRATION OF THE STRETCHING FORCE FROM ELON- 
GATION AND TRANSVERSE FLUCTUATIONS DSDNA MEASUREMENTS 

We would like to use partition function transformation properties under space rotation to 
derive the Strick et al. |5| formula, within the WLC model in the limit where the molecule 
contour L is larger than the persistence length A. This appendix is an updating of an 



unpublished internal note of the author 
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The thermal averaging over the various configurations of the molecular chain is performed 
here within the continuous version of the WLC model, using a technique inspired by quantum 
mechanics. It leads to the following expression for the partition function: 

Z(t , t 1; F, L) = £A(t • z)*„(ti • z) exp - (e n (a)^) , (A.l) 



A, 

where to, ti are the unit tangent vectors along the chain at the two ends so = and si = L, 
e(a) and \& n (t -z) are respectively eigen- values and eigen-functions of the WLC Hamiltonian 
Hwlc 0|- The partition function written above is clearly invariant upon simultaneous 
rotations of to,ti and F but it is modified by rotating F while keeping to,ti fixed or vice 
versa. If the experiments are performed under the condition L> A it is legitimate to make 
two simplifications in the right hand side of Eq. flA.ll) : 

i) The sum over n is dominated by the ground state contribution (n = 0), since the 
excited state contribution is strongly suppressed by a factor of the order of exp —(AeiL/A), 
where Ae x = e x — e ~ 1. 

ii) Keeping only the ground state term, the logarithm of the partition function reads as 
follows: 

lnZ(t ,ti,F,L) = —eo(a) (l (A.2) 

The term 0(j?) comes from the pref actor involving logarithms of the ground state wave 
function. 
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So we can conclude that the free energy of the molecule stays invariant when one rotates 
the force F, while maintaining fixed the two end tangent vectors to, ti, provided finite size 
effects of the order of ^ are ignored . 

We are now going to exploit this result by using the path integral form of the partition 
function: 

Z(t„, h, F, L) = J Vit) exp - ( ^^J^** ) , (A.3) 

where T>(t) is the path integral measure and Ebend is the elastic energy describing the re- 
sistance against bending. The stretching energy E stre tch is the potential energy associated 
with the uniform force F applied to the free end of the DNA chain: 

Estretch = -r(L) ■ F = - [ L ds t(s) ■ F , (A.4) 







where r(L) is the coordinate of the end point. Instead of taking as usual F along the z-axis, 
let us apply to F a rotation 1Z y (S) of angle 5 about the y-axis, leaving unchanged to and t%. 
The stretching energy reads, then , as follows: 

E st retch(5) = ~(z(L) cos 5 + x(L) sin 5)F. (A.5) 

By writing that the Taylor expansion of In Z(i , t 1; F, L) in the vicinity of 5 = vanishes 
term by term, we shall obtain a set of linear relations involving the moments of x(L) and 
y(L). Each moment is weighted by powers of / = such that the relations involve 
dimensionless quantities. In fact, the relevant relation is obtained by writing: 

.. 1 d 2 Z{5) ( E stretch {5)\ d 2 ( E stretch (5)\. 

= f 2 (x 2 (L))-f(z(L)), (A.6) 

where the r.h.s. has to be understood as a thermal average. If one uses the rotation 
invariant partition function given by eq. (jA.2[) . the l.h.s of the above equation vanishes, up 
to corrections of the order of j-. One arrives then immediately at the formula of Strick and 
5j which gives the stretching force F in term of the DNA elongation and the transverse 



al. 



fluctuations of the free end: 



(x 2 (L)) V L 
Note that the above derivation does not use any small angle approximation 



F= ^ L l k f + (A. 7) 



24 



We shall try to estimate the correcting terms under conditions such that -? is no longer 
negligible, say ^ ~ 0.2, while exp (—4) ~ 0.007 is still very small. The simplification i) is 
still legitimate, so we can use for the partition function the following approximate form: 

Z(t , ti, F, L) ~ * (to • z)* (ti • z) exp (-e {a)^) . (A.8) 



A 

If we apply the rotation 7Z y (5) to the force, the partition function of eq. (1A.8|) is then 
oc (cos(6'o — 5)) (cos(^i — 5)) . The computation of the finite size correcting term: 
A/. s (0q, 9\) = lim^o zjsj 9 ds*^ * s straightforward, but the final result is rather lengthy. The 
calibration formula corrected for finite size effects, valid for any value of 9q and 9\, is then 
readily obtained from eq. (1A.6j) : 

^w( i+4 " M) (4 + <»)' (A - 9) 

With the above writing, the finite size correction is proportional to A/((z(L)) a). We are 
going now to estimate the prefactor Af, s (9o,9i). It is of interest to quote the result in the 
simple case 9\ = 9 = : 



Let us first consider "large" force values : a » 1, remembering that a = 4 corresponds 
to F = 0.3 pN ! An approximate ground state wave function can be easily derived, together 
with the corresponding eigen-energy: 

\I/o(cos6>) oc exp — ^{yfa9 2 ) and e («) = — a + \fa . 

It leads immediately to the relative extension: (z(L)) /L ~ 1 — ( this value is very 
close to the exact one when a > 4). The correcting term Ay. s (# , #i) can be computed easily 
for arbitrary 9i angles: 

A/. S (0 , 0i) = -2v^ + a (9 + 9 ± ) 2 . (A.ll) 

The above formula shows clearly that something has to be said about the angles 9q and 9\. 
The angle #0 gives the direction of the initial strand, having a length of about 0.1 A, sticking 
out from the anchoring plate. The angle is partly determined, first, by the biological gluing 
process, second, by the fluctuating tension force induced by the rest of the chain. If the latter 
mechanism were the dominant one, a thermal average would have to be performed, using 
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the probability distribution oc \I/q(cos0o)- A similar analysis holds for the terminal strand, 
sticking into the magnetic bead. Computing the thermal average ({{Oo + #i) 2 ) = 2((6> ) 2 ) = 
2y / a, one gets A/. s (6'o,^i) = 0. ( This result turns out to be valid for arbitrary forces. It 
hinges upon the fact that thermal averaging over 6q and Q\ and the second derivative 
are commuting operators. The latter statement is not totally trivial since the wave function 
^/q(cos6) is restricted to the finite interval < 9 < tt.) The above considerations lead to 
the following approximate bounds for the finite size correction: 

a»l=> ~2y^<F f 2(L)) -1<0. (A.12) 

Let us apply the above result to a situation considered in section II, involving a "short" 
molecular chain A/L = 1/12 together with a stretching force a = 4 — > F = 0.3 pN . The 
large a asymptotic formulae give already a fair approximation. One sees immediately that 
the finite size correction stays below the 10 % level, when a > 4. 

In the small a limit, the ground state wave function is obtained by a first order pertur- 
bation calculation: 

^ O (cos60 = -^(1 + acos8 + 0(a 2 )). 
V 2 

Performing some simple algebra, one gets: 

^y^- = - a and A/. s (0o, 0i) = -« (cos(0 o ) + cos(0i)) + 0(a 2 ). 
L 3 

In the present low force limit the initial and final angle, #o an d #i, are weakly affected by 
the tension of the chain. To take into account the spatial constraints associated with the 
anchoring plate and the magnetic bead we impose the restriction < 6^ < tt/2. We arrive 
in this way at the following bounds: 

a « 1 ^_^< F J^L_ 1 < . 

Let us take as a typical low force case a = 1/3 =^> F ~ 1/40 pN. Using the relation 
((cos9)) = (z(L))/L = 2/9, one gets for the average angle between the running tangent 
vector and the force direction a value of about 1.3 rad. If one takes A/L ~ 1/10, the 
correction can indeed reach values of the order of unity. However, if one deals instead with 
a "long" molecule with L = 300 A, as in the experiments analyzed in ref. 10J, one gets 
3A/(aL) = 3 x 10~ 2 . The formula flA.7[) is then rather accurate, despite the fact that it 
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was originally derived within a small angle approximation. 
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